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ABSTRACT 


This paper is an attempt to develop a method for 
determining the slamming pressures and loads on ship hulls, 
when the relative vertical velocity relationship between 
the wave slope and keel is specified. A basic two-dimensional 
theory for hydrodynamic impact is utilized which avoids 
linearization of the free surface. The mathematical model 
for the theoretical developement is a wedge. <A method of 
analysis is presented which yields the free surface shape 
and pressures associated with symmetrical penetration of a 
fluid surface by wedges of varying deadrise angle. fhe 
anaylsis is based on the principles of continuity of fluid 
flow and similarity. The principle of roman ah LY maker use 
of ‘the fact. .that the volume of displaced: water above the 
original water line. fhe prisciple of similarity accounts 
for the proportional expansion of the free surface with con- 
tinued penetration of the body. The free surface (consisting 
of both wave and spray) is considered to be increasing in 
such 4 manner that the slope and velocity of the surface is 
determined by the position and penetration velocity of the 
impacting body. Once the correct free surface shape has 
been determined by a graphical iterative method, for any 
instant of time, the associated velocity potential and 
pressure can be determined numerically for various points 
on the body. 


Since only V-shaped ship sections can be approximated by 
wedges, a method is developed by which the pressures on a 
U-shaped ship section of a low deadrise angle can be calculated. 
Results are presented which resulted from calculations 
carried out on a ship with low deadrise angles in the foward 
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sections. For comparison, a calculation was performed 

for a Liberty ship section which could be satisfactorily 
approximated by a straight wedge shape. Also, the integrated 
pressure load was determined for several instants of 

time for the ship being considered. 


The results of the calculations seemed to give 
realistic values for the pressure distribution. However, 
the numerical calculations, as developed in this paper, 
are quite tedious. Computer solutions need to be developed. 
Also, the approximation of ship sections of low deadrise 
angles by blunt nose wedges is questionable, since the 
spray root thickness has to be approximated. fTherefore, 
the free surface shape in the spray root area is an 
approximation. The theoretical results should be closely 
checked by experiments. 
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MOTATION 





Points on the free surtace 
c Nominal wetted half width of the impacting body (2. )} 
g Acceleration due to gravity (ft/gec. , 

b Wave height (double amplitude) (ft. ) 

H Normal draft (ft,) 

L Shiv length (ft. ) 

© Crigin of coordinate avstems 

o} Edge of epray; origin for s 

P Pressure at any point in flow fleld (1b/ in”) 

Distance frors erigin to point in fluid (ft. )} 


r Wave amplitude (ft, } 

R Amplitude of compiex number 

a, & Distance along free surface curve from ite origin 
at the solid boundary (ft, } 

t Time (sec. } 

Te Perlod of ancounter between waves und ship 

u Velocity component in x-direction (ft/ cee) 
Hesgultant veiocity at any point in fluid 
“Vv” + we (ft/sec) 

v Velocity component in y-direction (ft/gec) 

Vv Shin speed (ft/ aec) 

Va Water surface velocity relative ic modal botiorm 
{or flat plate); a funetion of y > c¢ (ft/sec) 

Ww Velocity component in z-direction (ft/sec) 

x, ¥, & Linear displacements of center of gravity of ship along 
fined axes 

x, Y, Z Coordinate system fixed in ship 

Zn Heave amplitude 
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Relative position of saip keel with respect to water 
surface (ht) 


Relative velocity of keal with respect to water surface; 


penetration velocity (ft/ zea) 
Angle between velocity vector, U, asd radius rv 
Deadrise angle 


Angle between radiai velocity and the surface velocity 


past the radius 

opray root thickness; paage angle for heave 
FPheas¢ angie for pltca 

2~coordinate of particular fluid particle 
Wave leagith (Ri 

y-coordinate of a particular fluid particle 


Mass density of Muid (eae. 


Pitch angle 
Pitch amplitude 


Fregueney of encounter between the ship and waves 
{rad/ sec) 
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1, DN PRODUCTION 


In recent years there has been a distinct trend toward increasing 
the speed and displacement of ships. ©uite often a ship's speed ia not 
governed by ita installed horsepower, but by its response to the sea- 
way. The term ‘response '' includes not only the ship's motion, but al- 
go the effect on the ship's structure caused by the motion. Cne possible 
result of a ship's motion in a seaway is the slamming force exerted on 
the forward portion of the hull structure. The force ia uaualiy caused 
by the emergence of the bow and subsequent impact with the water on 
reentry. The ship may vibrate for some time after such an irapact, 
and plates on the bottom shell may be damaged. If the impacts are 
severe enough, serious structural damage may occur, and even if the 
blows are of smail magnitude, but frequent, the hull will possibly be 
weakened by fatigue. 


The earliest theoretical work on slar:eming loads resulted from 
studies on the impacts sustained by seaplanesa when they landed, The 
quantitative evaluation of the impact forces was based on the “ expand-~- 
ing plate’ and ‘spray root’ theories. The expanding plate theory origi- 
nated with Von Karmen in 1929 (1). By thie theory, the water flow 
pattern around a wedge penetrating the water surface is taken to be equi- 
valent at each instant to the flow about a flat plate of the same width as 
the wetted width of the wedge. From this enalysis, the velocities, accel- 
erations, and resultant pressures and hydrodynamic reactions are ob- 
tained, Apparently Von Karman's theory is good for low deadrise angles, 
but is only a crude approximation at higher angles (2). By the apray root 
theory, developed by Wagner in 1932, the local flow phenomena occuring 
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at the edge of the wetted area is described. The high preseures that 
exist near the edge of the wetted area during impact are connected 
with the formation of a apray jet in which the water is quickly accel- 
erated to a high velocity. 


Seaplane impact theories heve been adapted to the problems of 
ship slamming (3, 4) in en attempt to predict the pressure response 
of a given hull form, Experiments have been devised and carried out 
(4, 5, 6) utilizing ship models. The experimentationz have heen carried 
out under conditions allowing three dimensional flow about the ship hull, 
which are obviously the cond 
dimensional flow analysis for e general curved body ig much more 
complex than that of two dimensional flow. Therefore, theoretical 
work has been primarily concerned with the flew akout two dimensional 
ship sections of finite length, The slamming loads for the individual 
sections have been combined by strip thecry to obtain the pressure 
distribution on the total area of the ship bottom affected by the slam. 





tions met by ships at ses. However, three 


There is 2 pressing need for the extension of the available two- 
dimensional water entry theories, In the expanding plate theory, the 
free surface is represented by a gero potential line, The impact of 
ship hull shape bodies on a water surface will usually result in the 
formation of a curved free surface, and a spray root. Theoretically 
the presence of a highly curved free surface may noi be represented 
by a atraight sero potential Line. It seems Hkely that the free surface 
ig nearly a horizontal line only for bodies of high deadrise angles, and 
that the spray root generated at low deadrise angles disputes the assump- 
tion of a horizontal sera potential line. In general, for any of the ‘fit- 
ting’ techniques (ellipse and eirele fitting ae well as the expanding 
plate), the agreement between theory and experiment nas been rather 
poor (2), even for simple wedge zhaped bodies. Therefore, @ more 
exact solution of the hydrodynamic impact problem would be desirable. 
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The object of this work is te formulate ea basic two-dimensional 
theory for hydrodynamic impact which ayoids linearization of the free- 
surface boundary condition, and whieh can be adapied te the ship im- 
pact problem. The theory ia based largely on the work of Pierson (7) 
who developed an iterative numerical goluticn for a wadge entering 
the water af constant velocity. While a wedge is uaed as the rmathe- 
matical model in this paper, the velocity is mot considered to be con- 








The problem of determing slamming pressures actuaily consists 
of two parte; first, the description of the ship motion in & seaway which 
leads to slamming, and second, the description of the pressures and 
forces on the mull caused by sleammiag. Ship motions are not treated 
here, but an expression for saip motion is needed to determines, at any 
instamt, the relative veriicai velocity of the ship with respect to the 
water surface, Such an expregsion was cbiained fram sxperimental 
work for a Liberty type ahip, discussed in reference (3). 


Because of the difficulty of analytically describing the flow about 
an arbitrary conves body, the expressions for the velocity potential de- 
riveé in the paper are based on the water entry of a wedge, The dead- 
rise angle of the wedge ia taken to be approximately the same as the 
deadrise angle of the ship section under consideration, While the above 
assumption aeppeara tc: be somewha! groas, in the ight of past expert- 
mental and theoretical work it does net seer: te he bed. Both experi- 
mental and theoretical calewintions have shown that the highest trang - 
verge preseure distributions occur om a given ehip section shortly after 
it hits the water (3, 6, 93. The treoraticai calculations were based on 
the expanding plate theory, with a modification for the mull section cur- 
vature and added mase of water, However, the portion of the hull tat 
experiences the largest pressure loading is the bottom portion whic, 
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for most ships, is reughiy emttivalent to a wedge of low denirige angie, 
or can at least be agproulmated (at least for the period just before and 
after slamming) by a wedge of « specified angle. It is uot sufficient 
just te compute the preasure distrloution at the instant of slamming, 
but pressure distributions for several lostania before end after slam - 
ming are aleo required for a complete picture of the pressure loading. 
From a structural viewpoint, it ie the pressure icading and ita time 
variation which are important, and not peak presaures which last only 
for a few micro-seconds. 





The development of the theory is presented in Chapter 0, The re- 
sulia of the theory, as applied to a given ship for which other theoretical 
tad experimental dats are available, are presented and discussed in 
Chapters If and IV. 
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i, FRCCEDURE 


A, General The theoretical calculations are carried cut for a ship 
approximately the sis¢ of a liberty ship, The actual ship Umea were 
calculated from the series of sectiona for which Professor F. M. Lewis 
obtained his values of inertia couffictents ($8). Professor J, £. Kerwin 
obtained the offsets and polynomial expressions for the ship lines by 

the use of the 1050 [BM computer. Also, a calculation is performed 

for comparison for a section of the Liberty shin used in references 

(3S, 4}. 


B. Determination of Ship Motier The coordinate system used for the 
determination of ship motion is the game one usetl ia reference (10), 
The coordinate origin is taken to ke at midships. The (z, ¥, u} coordi- 
nates represent linear displacements of the center of gravity slong 
fixed axes. In thia paper, only heaving and pitching motions are eon- 
sidered, therefore only the ¢ coordinate is cf interest. The z axis is 
vertical and 2 is positive upward, Of the three angular displacements 
(roll, pitch, and yaw), only the piteh angle is of interest. Fositive 
indicates bow up. Also, the (X, Y, 2) coordinate system, which (e re- 
lated to the ship's principal axes, is introduced. The 2-coordinate is 
positive up, the 4-coordinate positive forward, and the Y-coordinate 
positive to port. The pertinent coordinates and their positive directions 
are shown in Fig. 1. 


in order to obtain a complete picture of the ship motion at any 
instant of time, it is necessary to kuow the instantaneous draft and re~ 
lative velocity between the keel sud water surface for say station, 


The pitching motion is asaumed to be 


Y =. cas (Het + ¢ ) 
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PeGUre 1 PRINCIPAL COORDINATES AND AXES OF 
REFERENCE 


Undisturbed surface 


; 





Become RELATIVE ELEVATION OF KEEL ABOVE 
DISTURBED WATER SURFACE 








and the heave is asgumed to be 
#5 cos {Wot +d). 


The wave elevation at any station along the ship is given by 





20K Oi. 
T * Yom 298 6 r + Wet) 


where We = a (V + 3 ) = period of encounter between ship and waves 
in head seas. 


The relative position of the keel above the disturbed water surface 
is given by: 
2° -#- XY +H (i) 
for any station (see Figure II}. Emergence is associated with a, < 6. 


The relative velocity between the ship aad wave at a distance A 
from amidships is given by 


‘2 T-2-x (2) 


where? = SPs - wor sin (wet + “25 .), 





end (J, is the absolute angular frequency of the wave and equals ae 





e+ 8 we vin (Wet +»), 


i; ° we sin (Wet +¢). 


From reference (11), experimental values of 3 and were obtained 
for a 5.5 ft. model of a Liberty ship in regular waves of a specified 
height. The values given for the model are tne same for the full scale 
ship, with the exception of the heave amplitude, which ia scaled up for 
the ship by a ratio of the lengths. The pertinent variables for the mo- 
del, full geale ship, and equivalent Lewis ship are tabulaied in Table I. 


The calculations of the relative velocity between the keel aad the 
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water by Bq. (2) will result in a sinusoidal variation of velocity for 
any station. While this ia not correct, it is shown in reference (3) 
that thia is a good approximation for the relative velocity of a ship 
subject to slamming motion in regular waves. In other words, the 
change in relative velocity due to slamming loads is very small, 


Slamming, ag generally defined, ia the sudden change of the 
acceleration of the ship, While the relative velocity dees not deviate 
_ much from that predicted by (2), the influence of slamming on the accel- 
eration {a quite marked (3). However, in this study of clemrming, the 
pressures are governed by the velocity and penetration of the keel in- 
stead of the acceleration, For the ship hulle under consideration, the 
portion’ of the ship between stations 6 and 5 is subject to slamming. 
From reference (3) the instant of slamming is determined to be in the 
vicinity of t = 1.20 sec. The penetration, and relative velocities of ata- 
tions 0 through 5 are graphed for the interval from t = 1.30 sec. to 
t= 1,44 sec, im Figs. Hl and IV. These values will be used to deter- 
mine the slamming pressures on the forward portions of the ship. The 
values used to plot the curves in Figures II and IV were computed from 
equations (la) and (1b) which resulted when the amplitudes and phase 
angles from Table I were subetituted into equations (1) and (2). 
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a. 12.8 cos (Wet + 385 } - 3,35 cos (wat -0, §23)-0, 122Xcos(watr+d. 26 


(la) 
a * -9.8 winiugtt Soe )-8. 35 ein (ver - .628)- . 128 Xsin (cuet + . 252) 
(ib) 


C. Theoretical Analysis 

Conditions of the Probiem 
i. The problem is treated as a two-dimensional one, The velociiy po- 
tential and pressure distribution will be calculated for two-dimengional 
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Model ship type ship 
Renmei (BP). .L, ft. 416.0 420.0 
a a 
j 
Wave length a 379.0 382.0 


Fiteh phase , degrees 14.44 14,44 14,44 





Frequency of encounter Qk, rae 8.74 PaO: mere 
Oe a2 6.28 6.28 
6.36 OL is OP: 
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sections (each a station spacing in length), and then combined by strip 
theory to yield the pressure load on the forefoot of the shin. 


2) Bach impacting section is assumed to be perfectly rigid and ic be 
symmetrical about the vertical centerline plane. 


3) The velocity of penetration ia not constant, but varies according 
to equation (2). 


4) The flow field is assumed to be irrotational (a necessary condition 
for the existence of 2 velocity potential). 


S$) The fluid is assumed {o be frictionless and incempressible, and 
the velocity of penetration is assumed to be high enough for the effects 
of gravity and surface tension io be neglected, 


6) On the basis of the assiunption of irrotational flow, the particles 
of the eriginal surface remain on the surface and uo new particles are 
added to the surface (12}. Sime the fluid is incompressible, the arc 
length along the surface from. the impacting bedy te any given particle 
is constant and independent of time, 


The preceding conditions deo not define the shape of the free aur- 
face. However, for a wedge shaped body, waich ia the mathematical 
model used for the theory, two general conditions can be formulated 
which will define the shape of the free surfaces 


a) Continuity: The fluid displaced by the body must appear above 
the original water surface in the form of a wave and/or spray. The arc 
length of the surface must be constant. 





b) Similarity; Sinee the impacting body is symmetrical, and the 
immersion is normal to the surface, the field of flow is divided inte two 
haives by the centerline of symmetry. By neglecting the gravity forces, 
the surface shape and dynamic state of the field of flow can be defined 


entirely in terms of the penetration and velecity of the body (7). If the 
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penetration velocity ia constant, the flow fleid and surface shape must 
be geometrically and dynamically similar at alltimes, However, tn 
the present case, the penetration velocity is not constant, and there- 
fore the surface shape is assumed to he only geometrically similar at 
different times. 


CONTOVUITY The general shape of the free surface due to penetra- 
tion by a wedge shaped body ia shown in Fig. V. Ag shown in the figure, 
the free surface can be thought of as being made up of spray and a wave. 


The two merge together in the ‘spray root’ area, In determining the free 
surface shape, the wave rise (3) is given by equations 11 and 13 of refer- 


ence (13). 
ae LUE E sin? © -y) (3) 


This is the came expressica which gives the surface rise in the "expand- 
ing plate’ theory. The wetted half width (<) is the point on the body which 
is wetted by the piled up water in the absence of spray. The wetted half 
width at the original water level is ae (Wagner; see reference 14). 
The continuity condition ia satisfied (f Area A ~ Area B, and ifs « y. 
However, the continuity condition iteelf is uot sufficient to determine 
the free surface shape aince any thickness of sprav may be associated 
with an unlimited variety of waves, and still not vioiate the continuity 
requirement, Also, the free surface in the spray roct area, where the 
surface weve and spray merge together, could assume any one cf many 


shapes, 


Similarity The principle of similarity for the proportional {ncrease of 
the wave and spray with continued penetration is the second considera- 
tien. Not caly must the displaced fluid appear above the original water 
surface, but the wave and spray must be increasing in such @ manner 
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FIGURE V FREE SURFACE DUE TO PENETRATION 
OF A WEDGE 
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that the shape and velocities are everywhere determined by the position 
and penetration velocity of the entering body. By a geometrical con- 
struction and integration process proposed by Wagner (eee reference 7), 

it ia possible to reiaie the surface shape and velocities, wo that for any 
spray thicknees, one correapending wave and spray-root shape can be 
determined which meets the requirements of continuity and similarity. 
Mowever, Wagner's proposal wes lor a constant entry velucity, the follew- 
ing theory ia developed for a variable penetration velocity, 


In an ideal fluid where gravity, viscosity, and surface tension are 
neglested, the shape of the free surface is uniquely determined by the 
ahape of the impacting body, and the sige of the shapes is determined 
by the amount of penetration . 


Conaider the several stages of the penetration process of a wedge 
shaped body in Figure VI, u#ing an y-% coordinate system. The penetra~ 
tion of the body into the fluid causes the initially flat surface to rise as 
a wave, and the nesrest particles are deflected out along the wedge sur- 
face as spray, The resulting aurface shape expands as the wedge pene- 
trates deeper. The resulting surface shaper are o'; @, by, o',& be, ete, 
For the same reintive position is the fluid and om the surface, tie pres ~ 
sure and velocity depend wien the surface shape, and must bear the same 
relationship to the penetration velocity at any instant of time (7), The 
above statement apelies to the relative position in the fluid, and not to 
a particular fluid particle. The peaitiens are indicated by of, , o's, O's 
and a, & , %;, etc. If the penstration velecity waa constant for the bedy, 
the velocity of the fluid at corresponding points would be equal. That is 
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" Of course this statement ie only true for a certain definite time inter- 
val efter impact. As the body penetraies deeper, much of the free gur~ 
face shape will be dissipated ag spray, and as the penetration velecity 
decreases, the influence: of gravity will be felt, 
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U “, a = U ™ K z where Ke represents & proportionality facter 
between the velocity U at positions a; , a,, etc., and the penetration 
velocity, &. However, for varying velocity, UL, f u,. ‘vu ie’ but the 


velocity at corresponding points will still be proportional to the pene- 


trating velocity. 
4,” Ke e (ta), °. = Be t (t,) 
..* a (ts), 0, “Ks (t, ) (4) 


However, it is not unreasonable to asaume thet ‘. * . ak +. a 

® 
a a. ete. Therefore, the velocity at corresponding pointe will be 
Girectly proportional to the instaniangous penetrating velocity. 


The velocity at any point on the free surface, such as a, or b,, 
is determined from consideration of this condition of similarity. Con- 
sider point b, for example. A particle on the surface at point b, is 
instantaneously raeving along the radius with a velocity U rb and at the 
same time sliding past the radius toward the spray region with a velocity 
U,. (see Fig. VII). 


The vector sum of the two motions are +U a! yields the instan- 
taneous resultant velocity u, of the point in question. 


The radial expansion velocity, Ue the velocity of the surface 
past the radius, Ue and the resultant velocity, U, are each directly 
proportional to the penetration velocity, Z. Therefore, fer constant 
entry velocity, or for an average velocity taken over a smal! time in- 
terval, the velocity compszents (U,,. US at @ position are constant, and 
the vector sum may be obtained from the radial distance and the arc 


length. 


Consider the point b., on the curve in Figure VI. The radiai 
distance r is given by: 
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FIGURE VI REPRESENTATION OF SIMILARITY 
CONDITIONS 
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PIGURE VII VELOCITY DIAGRAM 
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The vector 3 is constructed tangent to the surface at b, so that 
ea equals the arc length 0'b.,; t 


a» UL, 


The vector sum ia then: 


it 
i++ [ue (5) 


The individual fluid particles do not fellow the radiating lines, 
but remain at a fixed arc length from @' . The mation of a particle is 
therefore a function of lis relative poultion en the free surface, and is 
given by the following integral tie 


is =. -[va {Ga} 


zi 
% = [wa (8d) 
0 


Equations (6) are referred tc the ¥ and & axes in Figure Vi. The motien 
in the y direction in a time interval t ia expressed by 5 - 5, © being 
the original position of the particle. The symbol » represents the vert- 
ical elevation of a particle (z-direction). Now tre velecity of a particie 
{v, w) is te be determined for any point in its path, so that Equations 

(8) can be solved. There does not seem to be a simple analytic solution, 
but a graphical solution can easily be carried cut. However, the solution 
can be simplified by changing the variable of integration from Ume to arc 
length. This is done in the following manner. If the penetration velocity 
was constant, the arc length and times at two correspencing points on 
the frea surface could be related by: 
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{reference (7) ). However, for varying velocity at time t, U ss By wet . 


where Eb ig the proportionality factor which relates the velocity of the 
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surface sliding past the radius 0b tc penetration velocity. However, in 


the time interval from t; te tg, we cau assume thai U5 is constant and 


b 
equal to its average value. Actording to referenca (15), aa overage valve 


for a function, f+), can be defined as follows: 
A 


A £ {t) = Stew ty see a 


, + 2 8 
Therefore, at time ie , thie averages value of U ab can ve &eesumed 
to be constant over the interval i, to t.. For any tima, t, , 
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By this reasoning, an expression for arc length in terms of time can 


be obtained for any time interval, <¢: 
S ts 
2) 


5 
where t « the earlier time, and 5 * corresponding are length on the sur- 
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face at time t2. 


Ste 
Also, de - a as 
Substitution for dt in Equations (6a) and (6b) yieids 
Do — a 
5 fs (reanen (Bee a 
i 7 
qo | ds (Tb) 
Do 


Equations (7a) and (7b) can be solved numerically as follows: 


Referring to Fig. V1, tt is seen that the penetration velocity has 
approximately a linear variation over the time interval under considera- 
tlon. Since the particle velocities at corresponding points are directly 
proportional to the penetration velocity, they can also be assumed to 
vary linearly. Using this assumption, the integral in Equation (5) can 
be expressed ag 
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B+t- fu, asus (5) 
where UL is the average value of UL fer the time interval under question. 
From this, the average value of the velocity components, in Equations 
(Ta) and (7b) can be found (vt and wt). Therefore, a step by step inte- 
gration of Equations (7a) and (7b) can be carried out as shown in Appen- 
dix A, Any trial free surface shape which satisfies continuity can be 
used for the determination of the surface velocities. The velocities 

are then used to determine a new surface shape. This procedure is 
carried out until there is a negligible difference between assumed and 


computed values of surface shape. 


The integral equations (7a) and (7b) relate the surface rise at a 
point to all the points more distant from the body. Therefore, the path 
of each particle on the surface foliows the same related path, and the 
conditions for similarity are satisfied. It would appear that the integral 
equations reduce the number of spray-wave combinations to one for each 


spray thickness (7). 


The spray root thickness itself was shown in reference (7) to be 
related to wetted width (c) and deadrise angle (8): 
é 


$eei> (8) 


Potential on Free Surface 


Referring to Figure Vii, the tangential velocity, 8 along the 


- 
z & 


surface at a position, bis given by - U. +U., cosy . For points along 





b ‘. rb 
the surface (at a specific depth) cos Y = Gs 
Therefore, 
a¢ er 
38 Uap ® ve cost = eb ee ‘ds (8) 
5 ———— ——— a 
1 <= for a particular surface should not be confused with = . */s for 
ds das 
& constant position, such as b. 
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Since U_, and U_. are proportional to é . and Z can be assumed equal 
to a moan vadae over the time interval under consideration, Us and U.. 
ean be considered approximately constant over tke same time interval. 
A differential equaticn can be obtained by multiplying Equation (4) by 
tds. (tUs~ as, ¢t U~ Ws 


~tdé « eds - rdr 


Integration yields 
(a* ~ x? ) 


tg 5 


+A (10) 
Since € + Oats reo, A=§, and gd = ———— (143) 
which ia an approsimate expression for the velocity potential on the free 
surface with varying penetration velocity. 


Potential and Velocity Distributions on Rody Surface 





Asguming that the free purface shape for « given Geadrige angle 
is unique, and is a result of potential flow, the velocity and potential at 
any point can be represented by 8 complex analytic function: 


Gee Univ, waht iy 


The Cauchy integral equation may be applied to the complex func- 
tion as described in reference (7). By Cauchy's integrel formula, a value 
of am analytic function in a region is determined by its values on the boun- 
dary. There is ne choice of ways in which the function can be defined 
at points away from the boundery omce the function ie defined ou tae boun- 
dary. Every change of valves of the function at interior points must be ac~ 
companied by a chenge of ita values on the boundary if the function is to 
remain analytic (18), Tha imegral formula for the point within the boun-~ 


dary ia: 
£ (2} 
f(a)* 3 [2 as (12) 
C > 
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gis the complex coordinste of tre position im the field. Since pointe on 
the svlid boundary are being investigeaied, the bounding contour circles 
only halfway around the point su thai 27: becomes ci, Aleo, Ne} » u + iv 
* (2) on the boundary, and equation (12) becomes: 


: 1. | Ute) 
0 (s,) . i {3 > (18) 
C o 


The expression can be seperated intyw tae real and imaginary parta by 
placing the origin of coordinstes at the point % under considerstion, 
and writing # in the polar forra: 


e« Rel? » git 8+ ie 614) 


da Tg cine + ido] (35) 


Substitution of the expressias for # and dz inco Rouation (15) ylelds 
| 4 
Uigp=U, -we S [e - iv [a (ina) + ido] (18) 
Rearranging jarma gives 
iv + ey = [4a + vdoO +i | uae - vd (ink) | > 
é 


eo that tie real and imaginary paris may be separgied: 


% : tf vee » Ud (ink) (17a) 
i | 
U. ef d vd (isk) (i vb) 


The construction for performing the numerical (ategration ef the 
abave equations is shown in Figure VII, The complex coordinates were 
oriented wo that the y~axie (lrasginery aig) t¢ along the wedge, positive 
toward the spray tip. The x-axis (resi agis) ic normel to the body and 
positive owt. The potential on the body surface may be obtained by inte- 
grating the velocity along the wedge from the spray tin te the apex, 


Pressure Diatribution 


The general equation for fluid pressure ia s potential flow of hi- 
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POLAR COORDINATE SYSTEM FOR THE 
DETERMINATION OF THE VELOCITY 
POTENTIAL ALONG THE BODY SURFACE 


oe 





compressibie fivid, ia used to deterinine the presture distribution 
oa the bodly (12) 


£ ant -¥ + F (¢} (18) 


The gravity force (©) and extrangwus impulees (F (t) ) are cero ainee 
gravity has been neglected, and the Muid waz assumed to be at rest at 
infinity, Therefore Equation (18) becomes 


Pa e -& (19) 


The resultant field velocity (J) can be determined on the free surface 
(Equation ($) ), and on the bedy (Equation (17) ). Now it remains to de- 
termine an expression for a , 


Assume that the depth of penetration at time t; is given by the aver- 
age penetration velocity (a), in the interval t © oto%*t,, multiplied 
by & . This ts the game as asauming constant velocity penetration in 
the specified time interval. At any position, (7+, 2)” in the fluid, the 
potential is therefore a linear function of time in the interval to t; (7). 


For any time interval t, to t, agsuming %... 
| t 
ér," Or, = (20) 
However, for a fixed point (y, 2) in the flow field, the shove expression 


ie not correct, If the peint (y, 2} cotnridme with the position Pe ‘.. at 
144. | ep ara oot 
t, , at time? = t+ 2 ¢ tie poettias y+ ab? z. ig @ siamese 2 * ; 


from (y, =). For emeli ct, the resultant veetar velocity can be sgeurmed 
ennstant in the region, 


ALt * % Ry, a) = @ 





ee mE SEPA ERIE PB EE De BSP NAAR ABR PNRS EAE EINE ON Oh FISD 
) See Figere IX, (7 x) refers to the distance to the eame relative point 

in the fluid at a Yot penetration Pe Therefore, as ifime increases, Y, 

increases, | 


* 23 - 


LA wee wT salesredals at beew a duds Ritbomrngaers 
(21) Yow wat ao 


ts) war Tr: a+ + 
eoets ewe wee | (1) F) eweivquet erencrtize base | 4) errs pile oT 
i feey of of Opomeres cow bintt oa) bow Derralyen aeee and esivesy 
mmoved (01) woltacp® evwhewe tT . clint! 

wn 2-84 
eeeWs 660) eal ne SEAlowereh oA nos (1) yieactm bieit tundiveny edT 
~s2 @ salary OY wot .| (TI) oolieupll) yYrow eds no bas .( (3) mati) 
, -- tal Goarrge ae walarres 

cen alt a 6 5 sD ts etareny te Sp a ta ae 

vol iim , 0 * 1 ato + | Sewrrenat edt ai ALA) Ulonlew onthowinaeg one 
af aanaRnany Cieaior tanttems aulmeses Of couse of ai Gift . awd 
oe! Sait ew ai” Ls ct) cahiweg Yoo 2A Lovecte emit Retsegs ef 
7) of oF Lewontin. 8F G) auadd Do wncisuen? camel) & syotered! «) Liltmsiog 


_ b etownes Ut A Greet ont oe WT 
oie Rp 
(oa) a ate* 
A ew a Alieti wal) eal cl fe 4) cieg bee e et <evewol!l 
ae a” blown oth dike eabinaies (2 (0) Pe tam att eerie tee of 
4 omar 2 4) 2 Cen egt ea eet ata pee. 


came olen eendee etey Ene ne 4 Dae wt 2 eet 
.Dalgers ae +) fretenns 


ye ee 






















aay orunes use tlw anna edt 06 woes 





vey T eek 


ee 
ea eee A ros emary (6 40) i weet oot © 
See Ber 






Zs 
O f 
1h See, 
D= Dr 
Zo, 
af ay ane 
O —F6, (Y,2) yf 
K ORIGINAL WATER 
SURFACE 
Ze, a = 
Zo, =| EN a AE 9! 
f / 
£2 Las Gee 
€ 
z p> On, ts 
a (y, 2) 2a 





ORIGIVMAL WA&TER 
SURFACE 


£07 fp “ 
fo 6 == cy 
fp Le = PE Zx br 
us LS bs : 


FIGURE IX FIXED AND RELATIVE POSITION POINTS 
imine FLULD 


~26< 





U eos a is the velocity component elong . Therefore, 





Le ay, z) t, At ~ Oy, &) o @ & ‘ ~ UALS cosa - rt, (21) 


From Equation (26) 
= Brot. « at 
7 ~, | re, 4 ty “Tt, 
Equation (21) then becomesr 
a: a at $ on 
age b "Tes UAt ees a (32) 
Taking the limit of Squation (22) as At approaches 6: 
>. lim 48 lim 
st  Ate@ “At” Ae Te at +e? 
sy «(Uewee (23) 


BSquation (23) now applies to all points in the fluid. The pressure distribution 
ia given by subutituting Equation (33) in Equation (18). 


es »B ad. | 4 
e ae U cos a F (24) 


On the bedy surface, the pressure may be found by using the values of — , 
%, U, ande already determiaed. The total load on the body is obteined 
from the integration of the pressure across the projected wetted width, 


Application te Ship Sections 


There are two possible methods of solving the water entry provlem 
of bodies eof geserai shape by considering the fere surface shape. The 
first {s ta construct the free surface based om the boundary of the actual 
ship section aa it enterr the water. However, thie method is questionable 
since the similarity principle hax act been oreved te hold for a body of 
genera: shape. The second method is to aporoximate the ship section by 
a wedge-lke shape of an equivalent duadrise angle. 


The period of interest in a glamming investigation, for a certain 


~ BF ~ 


wT . geeks ceoqeeee Wincier On) Bi & 260 U 





CR) femoral: oie eM bn, gO th eee » 
(0) aovep® ” = 
ri Pe gt = BI gt rt 
i) # Boe f00 - eee t. 
Te 
{asad e- nae rs a. at 





we amet tte 
ee eee oD bk ee ol eateg Lin of ewdiyys won (C2) aatteupe 
Hi | esoTawyS 2 (EL) aattaup® peduiiiedwe qv sevig ai 
ea) Pane zrt--4 
Ye Ramer eal? qaeeer ed Aen? oct Quer Wuweeees od cosbure ybos adi aD 
eniihe a) (het ats pe Gaol late aT Sealargere® wewcla chne UY 











Geico (re came emt gatvioe te eimai: silascq ow? rus eve) 

ot! s@hee suAdiee owt Od parichranic qa pqemte Lesenes be epiiond Tes 
lecror vet Yo cpabeued wh fs Seed teeters reel wel foutinzos of 4! aertit 
Vita het G 6] Rete C1 Lee yO vee est Pwhee 1] as onlicwe gida 


40 Goltwee eyialh @c) (Melee af 4 Goothe Gtueea dT .equde lereoey 
Mae SAMaew ieelev lps ut by eqede edil-mybrw 0 


(tadtes ete) pocleglowrrn. god mowede @ 4) @weeiz! Yo bebteq aT 


Ship section, is the time pericd immediately after impact with the 
water, The maximum impact forces is produced elther at tre ingtant 

of impact, or immediately after entry. For the Liberty ship investi- 
gated in reference (3}, it wae found thet the slamming presanures were 
of little significance by the time the draft at station & waa 9,3 ft, In 
fact, the times interval cof significant pressures wae lees than 0.17 sec, 
For this problem, the time interval wuler investigation is 0,24 aec. 


For & hull section with sharp curvature, the zpray will not neces- 
sarily adhere to the boundary. From photogrepis of wedge impacts in 
referances (14, 12%), the spray, after it leaves the wedge curface, shoots 
off at an angle roughly equivalent to the deadrise angle {eee Fig. %). 
Therefore, it seems possible that V-chaped sections, or sections with 
fine lines, could be approximated by a wedge of a certain deadrigs angle. 
However, slamming is mere dangerous for sections wiih srnali deadrise 
angles. li la obvious that these sections could not aatisiactorily be re- 
presented by @ V-shaped wedge. 

For U-shaped sections with emall deadrise angles, the method 
for constructing the free surface ahape will nave to be modified. It is 
obvious that the wave rise cannot be estimated correctly by Equation (3). 
Figure XI is a sketch of a typical U-shaped aection, The pressure diatri- 
bution for such « section {se satisfactorily predicted by the expanding plate 
theory, except in the apray root region, from y= btey«c. Kis there~ 
fore reasonable to expect thet the pressure on the flat portica of the hull 
can he given by assuming the pressure distribution for a flat plate fromm 
y*Gtoy*b. The flow past this flat plate” can be used to estimate the 
wave height at the free surface, From (18) the distribution of the normal 


velocity fer a flat piste of width ga, at — y > a, is given by 
Ve 
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Vi, will be the velecity of the free surface particles, The weve height 
at & point y can be found by; 


t Ce 
, Cy) + vgs ( rioge (37) 
0 O y* 


For the wave height determizaticas, the dimension a ia taken te be 
equal to c', the width of the section immersed at the original water 
line (see Fig. X1). When constructing the free surface, the apray is 
aseumed to shoot off at an angle determined by the alope of the hull in 
the spray reot area. This is the game as assuming that the side of the 
Ship is straight, and that the impacting bedy is reughly equivalent tu a 
dlunt nosed wadge. Acsuming that the similarity principle still holds, 
the free surface is constructed by using the arc length and continuity 
ronditions, and checking by Equation (8). 

For a U-shaped 2ection, tie velocity potential for the Mat portion 
of the hull (from v « @ to y = b) can be asaumed to be that for irrotational 
flow about a flet plate, From (19) the complex velocity potential past a 
plate perpendicular to the flow {a given by: 


¢ es +-iv/y-s* «@ -,ya'- 3 (23) 
ation (17) can be used to compute the velocity potential on the hall 


from y » » to y > c, The two velocity potentials sxould agree at y = b 
for any specific time. (A sample caleulation is discugsed in Appendix A.) 





The calculations, from which the results {n thie paper were de- 
termined, were carried out for the Lewis type ship diacussed previously 
(se Figures XJ). 
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I, RESULTS 


Included in this section ere the curves of the theoretical vressure 
distribution for the chip under consideration. These curves are shown 
in Fi gures XLI- AV, Alee, one transverse section of the Liberty ship, 
used as a model in reference (3), was imvestigated. The pressure distri- 
bution for this section is shown in Figure XVII. This gection waa quite 
V-shaped and had a deadrige angle cf 2 420°. Straight wedge theory was 
used for calculation of the pressure distribution, end the results are com- 
pared with the results given in reference (3) for the sama section, 


A. sample calculation for station 3 of the Lewis section ship ip die- 
eussed in Appendix A, The tangential velocity iv? om the Mat partion of 
the hull was calculated by two methods {Teble IV}, The firs? method was 
by the assumption of flow over a flat plate. This is the method used for 
all of the low deadrise sections under consideration. However, for com- 
parison, the tangential velocity on the bottom of the hull was calculated 
by use of the Cauchy integral formula and equatica (17). The results 
cheek closely with those frem the flat plate Tlow assumption. 





The pressures load on the forefoot was determined by strip theory 
for scverai instants of time. The integrated lead distribution is shown 


in Figure AVI. 
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IV_DISCUSSION OF RESULTS 


A, Introduction 


Before discussing the results, the theoretical assumptions made 
will be restated, 


1) Cnly two dimensional sections were considered. 

2) Each impacting section was assumed to be perfectly rigid, and to 

be symmetrical about the centerline plane. 

3) While the velocity of penetration was net constant, for the small time 
interval under consideration, it was considered to vary linearly and aver- 
age values were used. 

4) The flow field is assumed to be irrotational. 

5) The fluid is considered to be frictionless and incompressible, and 

the effects of gravity and surface tension are neglected. 

6) On the basis of irrotational flow, the particles of the original surface 
are considered to remain in the free surface. The surface deflected by 
the impacting body ia assumed to have a constant arc length. 

7) The theoretical procedure was based on a mathematical model of a 
wedge with constant entry velocity. The theory was applied to the water 
entry of actual ship sections by assuming thai the impacting sections 
were approximately blunt nose wedges. That is, the free syrface formed 
during the initial impact stages was considered to be equivalent to that 
formed by a wedge of an equivalent deadrise angle. The spray was con- 
sidered to shoot off at an angle determined by the slope of the hull at the | 
wetted width. For the hull sections where the deadrise angle was approxi- 
mately zero, the velocity potential was determined by the expression for 


flow along a flat plate. 
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8) The ship motion expressions were determined experimentally for a 
Liberty ship of approximately the same dimensions as the ship found 
from Lewis sections. 


B. Determination of Pressure Distributions 


The ship motion expression for the relative velocity did not include 
the effect of ship forward speed. For a small pitch angle, the total re- 
lative normal velocity of the keel is: 


Vv, UY? + KY +2-¢r 


where U is the effect of the forward speed on the normal velocity. 

If the pitch angle is smail at the instant of slamming, U can be safely 
neglected, Att = 1,284 sec., the first time for calculations, the pitch 
angle was approximately 0.2° . Therefore, the contribution of U to the 


relative velocity was not important. 


Figures XIII to A V represent the pressure distribution on the fore- 
foot of the ship at approximately the time of slamming. The pressures 
were computed for all stations from 1 to 5 which are immersed at the 
particular time. The pressure distributions are not really sufficient 
to give ge definite value to the maxirnum pressure to which the ship is 
subjected. 


The pressure distribution curves are similar to the curves ob- 
tained in reference (3) for the Liberty ship. While the peaks on the 
curves obtained in this paper were lower, this is probably not significant. 
The peak pressures are of extremely short duration, and are relatively 
unimportant, For example, at t * 1,270 station 3 was out of the water, 
but .0164 sec. later it had a peak pressure of 550 psi at y = &.3' , and 
.016 sec. later there was a pressure of approximately 130 pei at y = 5.38’. 
However, it does seem that the average presgure distribution at a station 


varies directly as the peak pressure. 
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An attempt was made to determine the pressure at station 3 at 
a slightly earlier time, but the construction of the free surface was 
extremely difficult to perform, and the author had no reagon to believe 
the results were correct. It is quite possible, however, that the peak 
pressure might be even higher at station 3 at an earlier instant of time. 


Referring to Figure AVII the pressures on station 2 of the Liberty 
ship were uniformly lower than those determined in reference (3) by 
the expanding plate theory. This result is probably the more accurate, 
since the ship section could accurately be represented by a wedge during 
the time of initial entry into the water, The wedge assumption would not 
have been valid at a draft of 6 f., or greater, but by this time the im- 
pact preseures have been greatly reduced, From the results given in 
reference (7), it is to be expected that the computation of pressure 
distributions on wedge-shaped bodies, by considering the free surface, 
will give somewhat lower pressures than those calculated by the ex- 





anding plate theory. 


The integrated glam:ming load distribution in tona/ft, of length, 
is shown in Figure XVI. At best, thease curves are crude approximations, 
since each curve was determined by oaly 3 or 4 points. A much better 
method would be to compute the pressure distribution at é and i sta- 
tions in the area of slamming. This was not done in this paper since 
the offsets for fractional stations were not available. However, even 
with the points available, a good idea can be obtained as to the total 


slamming load on the forefoot of the ship at the indicated times, 
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VY. CONCLUSIONS 


1) The method of analysis given in this paper is at least a reasonable 
approximation since the curves of pressure distribution are of the cor- 
rect shape. 


2) The procedure for determining the pressure distribution by con- 
sidering the free suriace ghape is a long and tedious process. There 
are inherent sources of error in the construction of the free surface, 
and in the graphical and numericai solutions for the pressures, even 
for a simple wedge. The procedure, when adapted to a ship form, is 
even more questionable. However, there does ssem to be some merit 
in the procedure, and it is concluded that it can be refined considerably. 
The biggest improvement would be to develop a computer solution for 


the procedure, 


3) The wave rise is probably accurately given for an impacting body 
by either equation (3) or (27), depending on the shape of the body, How- 
ever, it is evident that equation (8) cannot be used to determine the spray 
root thickness for a ship section. When the free surface was conatructed 
for the blunt nosed sections, the spray root thickness had to be estimated. 
This waa possibly a large source of error. 


a) The determination of the free surface shape in the spray root area 
is quite critical, This ie where the particle velocities are the highest, 

end even a small error in the shape results in a large error in directions 
and magnitude of particle velocities. However, even with the possible 
sources of error in the spray root region, the pressures actually obtained 
for this region were as realistic as those determined by expanding plate 


theory. 
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5) The procedure described in this paper is only valid during the 
initial {mpact, and a short time afterward. After the body has pene- 
trated to a sizeable depth, the effect of gravity, and spray dissipation, 
will render consideration of the free surface uselees, However, the 
nighest pressures will, of course, occur at the time of impact or irnme- 
diately afterward, when the added mass of the body is changing most 


rapidly. 
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VI. RECCMM ENDATIONS 
1) Experimental and additional theorsatical work sheuld be performed 
to determine if similarity principles can be applied tc ship sections of 
general ghape. Particular attention should be given to determining how 
the expansion of the free surface is determined by the shape of the im- 
pacting body. Ag a start, experiments could be devised which utilized 
impacting wedges of various deadrise angles aad a high speed camera. 
The spray root thickness and the angle ait which the epray leaves the 
body should be thoroughly investigated. 


3) Computer procedures should be developed fox the solution of 
hydrodynamic impact problems involving a non-sero velocity petential 
free surface, 


3} Accurate theoretical and experimental methods should be developed 
which will give the total slam toad acting on the forefoot of the ship at 
any instant of time. 


4) The effects on slamming loads of huil flexibility and water com- 
pressibiilty should be investigated both theoretically and experimentally. 
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APPENDIX A 


Sample Calculations U-Shaped Section (Lewia Ship, Station 3), 
t= 1,284 sec.) 
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The firat estimate of the surface shape is computed by using 
equations (27), and the continuity principle. From the plotting of the 
wave rise, the wetted width (c) is determined to be 5.95 ft. (See Figure 
A-I). The slope of the hull section et y = ¢ is determined to be approxi- 
mately 63°. This is the angle at which the spray {s assumed to shoot 
off, The spray tip, 0', is iocated by measuring along the wave and 
wedge from y = 10’ a diatance s = y = 10' . (At the distance y = 15', 
the wave slope is quite small, sc that s = y is a cloge approximation. ) 
A spray thickness of f/c = .1% wan assumed. Once the spray and wave 
have been merged together, the spray tip is lecated more exactly by 
following the curved surface (9' is located on the line drawa tangent 
to the hull at y « c, at an angle of § = 63°). The area of the displeced 
water above the origine!l water line is made to equal the area of the 
submerged body below the original water line. Then equations (6) are 
used to check the position of the water particles, The steps of the num- 
erical integration for the surface chape are shown in Table HW, The 
new particle positions were determined by the following anpraximate 
forma of equations (6): 
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After the free surface shape has been determined, the velocity 
potential can be determined on the body. For the flat portion of the 
huli, from y = 0 to y > 4.7, the potential ig determined by equation (28). 
In this computation, ais taken to be equal to the wetted width at the 
time uader consideration, but the equation ig only used for the flat 
pertion of the hull (from y = 9 te y = 4.7), 


For the portion of the huil (from y « 4.7 to y = 6, 05) the velocity 
potential is computed by integrating the velocity along the hull surface. 
The tangential velocity (vi) slong tne hull is determined from equation 
(17). Four points are chosen on the curved portion of the hull for the 
determination of the velocity (see Figure A-i}. The ly and x axis are 
constructed tangent and perpendicular to the hull at each point under 
consideration. The velocity iv is determined by equation (27) by num- 
evical integration around the free surface: 


73 =) vad - =‘. u d (in®) 
The steps in the integration are shown in Table Til for point 2(z, v, andu 
refer to the position and velocities of points on the free surface), On the 
body, ¢. ig equal to 2 com } at any point. Knowing the component velo- 
eitles, the resultant velocity (U) can be found at any point. The velo- 
cities and velocity potential determined at the points on the bull, shown 
in Figure A-I, are listed in Table IV, 


Wow the pressure can be determined for equation (24); 


See S .% 

. on Y= U cos a 3 
where is the distance froma the origin ic the point on the body under 
consideration, and« is the angie between and U. The results of the 
solution of equation (24) for station 3 a¢ i = 1.284 are listed in Tawie IV. 
The time, t, luvelved in equation (24) is the elapsed time from impact 
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of the body. For the case under consideration, t, {impact time) = 1.270, 
and the computations were carried out irom time ty * 1,284 sec, There- 
fore, t = ty - t « ,014 sec. in equation (34). 
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TABLE IT 
NUMERICAL INTEGRATION FOR SURFACE SHAPE 
Station 3 t= 1.284 z, = 0.3 i = 26.0 0 = 63° 


S 13; wt wt wtAs wt /\s } 
s? s+ 


start s* new 





Area of submerged body = 1.54 ft 


Area of displaced water = 1.60 ft 





alGe 


27.5 
Zone 


1.90 
4.20 
6.70 
8.30 
10.55 
1055 





TABLE III 


Numerical Integration of Equation (17) 





e010: 
v, (point 2) = =-(600) + (910) = 481 ft/sec 


UW = 2, 00863 = 11.8 ft/sec 


TABLE IV 


Velocity, Velocity Potential, end Pressure Distribution 
on Station 3 of Lewis Ship(t = 1.284 sec) 


Vy D Vo Vo Uo p 
rac (Flat Plate)| (Cauchy Integral) ( Pst) 
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